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Abstract 

We describe an algorithm for the Rosenbluth Monte Carlo enu- 
meration of clusters and lattice animals. The method may also be 
used to calculate associated properties such as moments or perimeter 
multiplicities of the clusters. The new scheme is an extension of the 
Rosenbluth method for growing polymer chains and is a simplification 
of a scheme reported earlier by one of the authors. The algorithm may 
be used to estimate the number of distinct lattice animals on any lat- 
tice topology. The method is validated against exact and Monte Carlo 
enumerations for clusters up to size 50, on a two dimensional square 
lattice and three dimensional simple cubic lattice. The method may be 
readily adapted to yield Boltzmann weighted averages over clusters. 

1 Introduction 



The enumeration of lattice animals is an important problem in a variety of 
physical problems including nucleation percolation and branched poly- 
mers 0. A lattice animal is a cluster of connected sites on a lattice with 

*Materials Research Institute, Sheffield Hallam University, Pond Street, Sheffield, SI 
IWB, UK 

I'CoUoids and Rheology Unit, ICI Wihon, PO Box 90, Wilton, Middlesbrough, Cleve- 
land, TS90 8JE, UK 



1 



given sjTiimetry and dimensionality and we seek to enumerate all distinct 
animals with a given number of sites. Exact enumeration has been carried 
out for small lattice animals using a variety of methods 0,^,0 but the meth- 
ods become computationally prohibitive for large animals. Many techniques 
have been used to enumerate larger lattice animals including various Monte 
Carlo growth schemes a constant fugacity Monte Carlo method 0, 

an incomplete enumeration method [|l^ and reaction limited cluster-cluster 
aggregation 0. 

In the following paper we describe an improvement of a method proposed 
by one of the authors ||TT| which was based on an extension of the scheme 



proposed by Rosenbluth and Rosenbluth [0 for enumerating self avoiding 
polymer chains. The central problem in using the Rosenbluth scheme for lat- 
tice animal enumeration is calculating the degeneracy of the clusters which 
are generated. In the method proposed by Care, the cluster growth was mod- 
ified in a way which forced the degeneracy to be A^! where is the number 
of sites occupied by the lattice animal. However the resulting algorithm was 
fairly complicated to implement. An alternative method of correcting for 



the degeneracy had been proposed by Pratt |T^. In this latter scheme the 
correcting weight is more complicated to determine and must be recalculated 
at each stage of the cluster growth if results are sought at each cluster size. 
However the Pratt scheme does not require any restriction on the growth of 
the cluster. 



In this paper we show that there are a class of Rosenbluth like algorithms 
which yield a degeneracy of A^ and which are straightforward to implement. 
The method provides an estimate of the number of lattice animals and can 
also yield estimates of any other desired properties of the animals such as 
their radius of gyration or perimeter multiplicities 0. We describe and 
justify the algorithm in Section ^ and present results to illustrate the use of 
the method in Section |^. Conclusions are given in Section ^ 



2 Algorithm 



Any algorithm, suitable for the purpose of the enumeration of lattice animals 
using the Rosenbluth Monte Carlo approach, must satisfy two important 
criteria. First of all it has to be ergodic. That is to say, the algorithm should 
have a non zero probability of sampling any given cluster shape. The second 
criteria relates to the degeneracy that is associated with each cluster and 
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requires this to be determinable. This degeneracy arises from the number 
of different ways that the same cluster shape can be constructed by the 
algorithm. While it is easy to devise methods of growing clusters that meet 
the first requirement, the second condition is more difficult to satisfy. For 
many simple algorithms the calculation of the degeneracy, for every cluster, 
can be a more complex problem than the original task of enumerating the 
number of lattice animals. 



In the original Rosenbluth Monte Carlo approach of Care 11 , this dif- 



ficulty was overcome by ensuring that the degeneracy for all clusters of size 
was the same and equal to A^!. However, to achieve this result the al- 
gorithm had to employ a somewhat elaborate procedure. This made the 
implementation of the method rather complicated, as well as limiting its 
possible extension to enumeration of other type of clusters. Here we shall 
consider an alternative algorithm, which while satisfying both of the above 
criteria, is considerably simpler than the algorithm proposed by Care. In 
Section |2.1| we describe the algorithm in its most basic form, before proving 



in Section 2^ that the ergodicity and the degeneracy requirements are both 
met. In Section |2.3| we demonstrate how the basic algorithm can be further 
refined to improve its efficiency. 



2.1 Basic Algorithm 

Having chosen a suitable lattice on which the clusters are to be grown (square 
and simple cubic lattices were used in this study for 2D and 3D systems, 
respectively), a probability p of acceptance and q = (1 — p) of rejecting 
sites is specified. Although in principle any value of p between and 1 can 
be selected, the efficiency of the sampling process is largely dependent on a 
careful choice of this value, as will be discussed later. In addition, an ordered 
list of all neighbours of a site on the lattice is made. For example, for a 2D 
square lattice this might read (right, down, left, up). While the order initially 
chosen is arbitrary, it is essential that this remains the same throughout a 
given run. In the basic algorithm, once chosen, the probability p remains 



fixed during the Monte Carlo sampling procedure. However in Section ^ 
the effect of relaxing this requirement is discussed. 

We construct an ensemble of Ne clusters and for each of these calculate 
a weight factor which we subsequently use to calculate weighted averages of 
various cluster properties. For a property O of the clusters, the weighted 
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average is defined as 

. Ne 

<0>;y=— ^l^A (1) 



Ne , 



The weight associated with cluster a with sites is defined to be Wa = 
l/{di\iPa) where is the normahsed probabihty of growing the cluster and 
djy is a degeneracy equal to the number of ways of growing a particular 
cluster shape. It can be shown |]TT| that the weighted average can be used 



to estimate the number, ctv, of lattice animals of size N and other properties 



?2 . 



such as the average radius of gyration R%- 

E[< 1 >w] = CN (2) 



E[< Rl >w] = ^ Rnu = cnRn (3) 

{u=l} 



During the growth of each cluster we maintain a record of the sites which 
have been occupied, the sites which have been rejected and a 'last-in-first-out 
stack' of sites which is maintained according to the rules described below. 
Each cluster is grown as follows 



(i) . Starting from an initial position, the neighbours of this site are exam- 

ined one at a time according to the list specified above. An adjacent 
site is accepted with a probability p or else is rejected. 

(ii) . If the adjacent site is rejected, a note of this is made and the next 

neighbour in the list is considered. 

(iii) . If on the other hand it is accepted, then this becomes the current site 

and its position is added to top of a stack, as well as to a list of accepted 
sites. The examination of the sites is now resumed for the neighbours 
of this newly accepted site. Once again this is done in the strict order 
which was agreed at the start of the algorithm. 

(iv) . Sites that have already been accepted or rejected are no longer available 

for examination. Thus, if such a site is encountered, it is ignored and 
the examination is moved on to the next eligible neighbour in the list. 

(v) . If at any stage the current site has no more neighbours left, that is 

all its adjacent sites are already accepted or rejected, then the current 
position is moved back by one to the previous location. This will be 
the position below the current one in the stack. The current position is 
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removed from the top of the stack, though not from the hst of accepted 
sites. 

(vi) . The algorithm stops for one of the following two reasons. If ever the 

number of accepted sites reaches N, then the algorithm is immediately 
terminated. In this cluster of size N is successfully produced. 

Note that unlike some of the other common cluster growth algorithms 
P] , it is not necessary here for every neighbour of the generated cluster 
to be rejected. Some of these might still be unexamined before the 
algorithm terminates. The second way in which the algorithm stops is 
when it fails to produce a cluster of size N. In this case, the number of 
accepted sites will be M < A^, with all the neighbours of these M sites 
already having been rejected, leaving no eligible sites left for further 
examination. From step |(v)| , it is clear that in cases such as this, the 
current position would have returned to the starting location. 

(vii) . The probability of producing a cluster of size N, in a manner involving 

r rejections, is simply Hence the weight, Wa, associated with 

the growth of the cluster is given by 

W^ = l/idNp'^'-''\l-pY) (4) 

where the degeneracy, d^, is shown below to be exactly A^. Failed 
attempts have a zero weight associated with them. However they must 
be included in the weighted average of equation (|I|). 

(viii) . During the growth of a cluster of size A^, we may also collect data for all 

the clusters of size M where M < N. It must be remembered that the 
weights for these smaller clusters must be calculated with a degeneracy 
of M. 

A specific example is helpful in demonstrating the algorithm. Figure ^ 
displays a successful attempt in forming a cluster of size A^ = 4, on a square 
lattice. The order in which the neighbours were examined was chosen to 
be right, down, left and up. Let us now consider various steps involved 
in construction of this cluster in detail. Beginning from the initial position 
labelled cell one, the adjacent site to the right of this position is examined. In 
this case the site is rejected and the current position remains on the cell one. 
Such rejected cells are indicated by the letter X. The next neighbour in the 
list is the one below, labelled cell two. As it happens this is accepted. Thus, 
the current position moves to this site and its position is added to the top 
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Figure 1: Sequence of accepted sites leading to a cluster of size = 4. The 
sites examined but rejected along the way are indicated by X. In our notation 
this sequence can be represented by {0,1,1,0,0,0,1} 



of the stack, ahead of the position of cell one. The process of examining the 
neighbours is resumed for sites adjacent to cell two. Once again, following the 
strict order in the list, the site labelled three to the right of current position 
is considered first. This is also accepted and as before is placed at the top of 
the stack. At this stage the stack contains the positions of cells three, two 
and one, in that order. The current position is now cell three. 

The site to the right of this, followed by the one below, are tested and 
both rejected in succession. Since both the neighbours to the left ( ie cell 
one) and the one above have already been considered, the current position 
has no more eligible neighbours left to test. Therefore, following step 



above, site three is removed from the stack. This leaves the position of cell 
two at the top of the stack, making this the current position again. The cell 
two has two neighbours, the adjacent sites below and to the left, which are 
still unexamined. Of these, according to our agreed list, the site below takes 
precedent, but as shown in Figure ^ this is rejected. Current position remains 
on the cell two and the neighbouring site (cell labelled four) to the left of this 
position is tested. As it happens this is accepted. A cluster of the desired 
size A^ = 4 is achieved, bringing this particular attempt to a successful end. 

For the subsequent discussion, it is useful to represent a sequence of ac- 
ceptance and rejections by a series of 1 and 0. Thus, for the case shown in 



6 



Figure |1| we have {0,1,1,0,0,0,1}. Note that at any stage throughout a series, 
the position of the current site and that of the neighbour to be examined, 
relative to the starting cell, are entirely specified by the decisions that have 
been made so far. In other words, given a sequence of one and zeros we can 
determine precisely the shape of the cluster that was constructed. This is 
only possible because of the manner in which the neighbours of the current 
position are always tested in a strict pre-defined order. For an algorithm 
that considers the neighbouring sites at random, the same will clearly not be 
true. 

The procedure described above needs to be repeated a large number of 
times, to obtain the weights for the ensemble average defined in equation (p. 
In particular, using equation ^, the number of lattice animals of size can 
now be determined. 



2.2 Ergodicity and degeneracy of the algorithm 

Let us now discuss the issue of the ergodicity of the algorithm. We wish to 
see whether, starting from any particular site on a given cluster, a series of 
acceptance and rejections (1 and 0) can always be determined which leads to 
that cluster shape. We stress that we are not concerned about how probable 
such a sequence is likely to be, but merely that it exists. We can attempt 
to construct such a sequence by following the same rules as our algorithm 
described above, with one exception; we accept and reject each examined 
site according to whether it forms part of the target cluster shape or not. 
Obviously, in the original algorithm, each such move has a non zero chance 
of occurring, provided p is not set to zero or one. Since we only accept 
sites that belong to the cluster in question, it follows that if the sequence 
is successful then we would achieve the desired cluster shape. However, we 
might argue that for some choice of target cluster and starting position, a 
series started in this manner will always terminate prematurely. That is to 
say, it will inevitably lead to a failure, with only part of the required cluster 
having been constructed. Now, it is easy to see that this cannot be true. If 
the series fails, it implies that all the neighbouring sites of the sub-cluster 
formed so far are rejected. However, the rest of the cluster must be connected 
to this sub-cluster at some point. Hence, at very least, one neighbouring site 
of the sub-cluster must be part of the full cluster and could not have been 
rejected. Starting from any of the sites belonging to a cluster then, it is 
always possible to write down a sequence of one and zeros that will result 
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in the formation of that cluster. Similarly, considering every starting point 
on a cluster of size N, another implication of the above result is that the 
corresponding cluster shape can be generated in a minimum of at least N 
distinct ways. 

Next, we shall show that the degeneracy of a cluster of size in our 
algorithm is in fact exactly (unlike the original algorithm of Care |Tl[ 



which has a degeneracy of N\). Let us suppose that starting from a particular 
site on a given target cluster shape, our algorithm has two distinct ways of 
forming this cluster. Associated with each of these, a series of one and zeros 
can be written down, in the same manner as that indicated above. The 
two ways of constructing the cluster must necessarily begin to differ from 
each other at some stage along the sequence, where we will have a 1 in one 
case and a in the other. Now since up to this point the two series are 
identical, the site being examined at this stage will be the same for both 
cases. This is rejected in one sequence (hence 0) whereas it is accepted in 
the other (hence 1). It immediately follows that these two differing ways of 
constructing the cluster cannot result in the same shape. Using this result, 
together with previous one regarding the ergodicity of the algorithm, we are 
lead to conclude that, starting from a given site on a cluster, the algorithm 
has one and only one way of constructing the cluster. Hence, for a cluster of 
size A^, the degeneracy is simply A^. 



2.3 Refined algorithms 
2.3.1 Adjacent site stack 

During the growth of the cluster a stack can be constructed of all the sites 
which are adjacent to the cluster and still available for growth. When a new 
site is added to the cluster, its neighbours are inspected in the predetermined 
sequence and any available ones are added to the top of this stack. (Note 
that this stack differs from that discussed in Section ( p.l|) ). The choice of 
site to be occupied can be made from all the adjacent sites in a single Monte 
Carlo decision. Thus, if we consider the underlying process in the method 
described above, at each step there is a probability p of the site being accepted 
and a probability q = 1 — p of the site being rejected. We therefore need 
to generate a random number with the same distribution as the number 
of attempts needed to obtain an acceptance. The probability of making k 
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attempts of which only the last is successful, is 



Pk = V 



(5) 



where 1 < A; < oo and Yl'k'=i Pk = ^- In order to sample from this distribution 
we note that the associated cumulative distribution, Cm, is given by 



Cn 



q) = l-q' 



(6) 



k=l 



Hence if we generate a random number, 77, uniformly distributed in the range 
< 77 < 1, then a number m given by 



m = Int 



ln{ri) 
\n{q) 



+ 1 



(7) 



will have been drawn from the required distribution. Thus we generate the 
number m according to equation (|^) and use this to determine which site on 
the stack is selected, with m = 1 corresponding to the site at the top of the 
stack. If m > Nadj, where Nadj is the number of available adjacent sites, the 
cluster growth is terminated as explained in step (|(vi)| ) in Section ^]T[ All 
the adjacent sites lying above the chosen site in the stack are transferred into 
the list of rejected sites. The list of adjacent sites is then adjusted to include 
the new available sites adjacent to the recently accepted site. As before, it 
is crucial that these are added to the top of the list in the strict predefined 
order. 



2.3.2 Variable probability 



An apparent disadvantage of the methods so far described is that with fixed 
choice of probability, p, occasions arise when a cluster growth will terminate 
before reaching a cluster of size A^, simply because the Monte Carlo choice 
rejected all the neighbouring sites. This problem can be overcome if the 
value of p is allowed to vary as the cluster grows. The simplest method is to 
determine the number, Nadj, of available adjacent sites at each point in the 
cluster growth and select one of these sites with uniform probability. This 
effectively makes p = 1/ Nadj and thereby increases the chances of growing a 
cluster of size A^. Note that it is still possible for a cluster growth to become 
blocked. This happens when the chosen site is the one at the bottom of 
the current eligible neighbours list, thus causing all the other neighbouring 
sites in the list to be rejected in one step. If the newly accepted site has 



9 



itself no unexamined neighbours to add to the hst, the algorithm terminated 
prematurely. Modified in the manner described above the weight associated 
with a cluster is now 



TTAf T\Ti 



N 

rather than the expression given in equation @). 

However, when this variable probability method was tested it was found 
that although it reduced the number of rejected clusters, it was inefficient 
at sampling the space of possible clusters when compared with method de- 
scribed in section (|2.3.1|) . This inefficiency was measured by comparison of 



the standard deviation in the estimated cluster number for any given number 
of clusters in the sampling ensemble. It is thought that the inefficiency of 
the variable probability method arises because it gives too much weight to 
sites lower in the stack, yielding many non-representative clusters. It is pos- 
sible that this problem could be overcome by using a non-uniform sampling 
distribution (c/ |jll|) but this was not tested in this work and the method 
described in ( p.3.1|) was used to obtain the results described in Section (|]) . 



3 Results 



In order to test the algorithm described in Section (0) it was used to esti- 
mate the number of lattice animals on a square 2D lattice and a simple cubic 
3D lattice for which exact results are known up to certain sizes 0. Before 
collecting data it was necessary to determine the optimum value of the prob- 
ability p with which an adjacent site is accepted during the cluster growth. 
The effect of changing p on the estimated error in the number of clusters of 
size 50 on the 2D and 3D lattices can be seen in Figure 0. It can be seen that 
there is a fairly broad range of values of p for which the error is a minimum 
and a value oip = 0.6 was used to obtain the results described below for the 
2D lattice and 0.72 for the 3D lattice. The distribution of weights is log nor- 
mal and becomes highly skewed for large cluster sizes; this is a standard 
problem with Rosenbluth methods [Q. The minimum in the error achieved 
by the choice of the value of the probability p has the effect of minimising 
the variance of the distribution of the weights, W^. 

In Table |I| we present results obtained using the algorithm defined in 
section ^ using the adjacent site stack method of section ^]3| to enumerate 



clusters on a simple cubic 3D lattice for clusters up to size 50. The results 
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Figure 2: Percentage errors for clusters of size 50 

were obtained from an ensemble of 2.5 x 10^ clusters. The data took 3.3 
hours to collect on a R5000 Silicon Graphics workstation using code written 
in the language C but with no attempt to optimise the code. Only 30% of the 
clusters achieved a size of 50. The results are quoted together with a standard 
error, e*^**, calculated by breaking the data into 50 blocks and determining 
the variance of the block means for each cluster size. If the number of samples 
in each block is sufficient, it follows from the central limit theorem that the 
sampling distribution of the means should become reasonably symmetrical. 
We therefore also quote a skewness, ^, defined by |T5| 



/ 3/2 



(9) 



where rrii is the i moment about the mean of the sampling distribution. 
It is expected that ^ ~ 0.5 for a symmetrical distribution and ^ > 1 for a 
highly skew distribution. The statistic ^ should be treated with some caution 
since it is likely to be subject to considerable error because it involves the 
calculation of a third moment from a limited number of data points. 



Exact results are known for clusters up to size 13 
quote the values for the quantity x defined by 
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Figure 3: Variation of percentage error with cluster size. 



where cm is the number of clusters of size M and it can be seen that all the 
values of x ^ire 0(1). Hence we assume that e^'^* is an acceptable method 
of estimating the error in the method. However it is likely that the e^'^* will 
underestimate the true error if the distribution becomes more skew. We also 
quote in Table |I| the values of cat calculated by Lam |Q using a Monte Carlo 
incomplete enumeration method together with the error estimates reported 
for this method. 



In Table ^ we quote data collected from a square two dimensional lattice 
by collecting data from 2.5 x 10^ clusters up to size 50. This data only 
took 1.45 hours to collect but only 2% of the clusters achieved a size of 50. 
Comparison is given with exact results p up to clusters of size 19. The 
rate of growth of errors for the two and three dimensional data is shown 
in Figure ^ and it can be seen that the errors associated with the method 
diverge are beginning to diverge quite rapidly above clusters of size 50. This 
behaviour is to be expected with a technique which is based on sampling from 
a log normal distribution. In the previous paper equivalent results were 
obtained for clusters up to size 30 with approximately the same sample size. 
The improvement up to clusters of size 50 obtained by the new method arises 
because the weight associated with clusters of a certain size is generated from 
roughly half as many random numbers. This effectively halves the standard 
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deviation of the log normal distribution of the weights and allows larger 
clusters to be sampled before the method becomes unusable. 



4 Conclusions 



We have described a simple Rosenbluth algorithm for the Monte Carlo enu- 
meration of lattice animals and clusters which can be applied to any lattice 
topology. A merit of the scheme is that for thermal systems it may be easily 
adapted to include Boltzmann weightings following, for example, the argu- 



ments used by Siepmann at al [0 in the development of the configurational 
bias technique. Similarly, the method can be applied to calculation of the av- 
eraged properties of a cluster of a given size, in the site percolation problem. 
In this case we have 



< (1 - pyo >w _ e£i - PY'-Oa 
<{i-py>w Ea=iV^o.{i - py- 



<o>-- 



(11) 



where P is the probability of site occupation in the percolation problem of 
interest and ta the number of perimeter sites of the cluster a. Preliminary 
results also indicate that the method may be useful in the study of the 
adsorption of clusters onto solid surfaces. A possible numerical limitation of 
the method arises from the highly skew probability distribution of Rosenbluth 
weights which occurs for large cluster sizes. However the method presented 
in this work is able to work to considerably higher cluster sizes than the one 



described in [11] before this becomes a problem. 
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0.16 




1.61 




-0.11 


24 


1.882x10^9 




1.877x1019 


0.19 




1.68 




0.16 


25 


1.485x10^0 




1.480x1020 


0.21 




1.70 




-0.02 


26 


1.169x10^1 




1.168x1021 


0.21 




1.75 




-0.11 


27 


9.214x10^1 




9.209x1021 


0.20 




1.81 




0.06 


28 


7.316x10^2 




7.290x1022 


0.21 




1.88 




0.18 


29 


5.790x10^3 




5.786x1023 


0.24 




1.96 




-0.12 


30 


4.600x1024 




4.610x1024 


0.25 




2.01 




0.44 



Table 1: Table continued on next page 
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37 


9.631x10^0 






0.39 








1.08 


38 


7.691x10^1 






0.35 








0.18 


39 


6.203x10^2 






0.40 








0.27 


40 


4.984x10^3 






0.45 








0.54 


41 


3.999x10^4 






0.43 








0.35 


42 


3.205x10^5 






0.46 








0.23 


43 


2.605x10^6 






0.49 








0.35 


44 


2.100x10^^ 






0.62 








2.32 


45 


1.684x10^8 






0.71 








0.43 


46 


1.353x10^9 






0.69 








0.65 


47 


1.087x10^° 






0.58 








0.36 


48 


8.892x10^0 






0.68 








0.53 


49 


7.223x10^1 






0.79 








0.02 


50 


5.789x10^2 






0.75 








0.78 



Table 1: Continued:- Degenerate Rosenbluth estimate of the number of lat- 
tice animals of size iV on a three dimensional square lattice using 2.5 x 10'' 
sample clusters, each grown to iV = 50 with p = 0.72; exact values from 0; 
estimated values and associated errors from incomplete enumeration method 
of Lam 0; calculation of error estimate described in text; 'true' error is frac- 
tional difference of Rosenbluth estimate and exact value; x ^^^^ ^ ^ire defined 
in the text. 
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N 


Rosenbluth 
estimate 


Exact 
value 


% error 


True 
%error 


X 




2 


1.999x10° 


2 










3 


6.000x10° 


6 


0.02 


0.01 


0.22 


-0.48 


4 


1.900x10^ 


19 


0.03 


0.00 


0.00 


-0.65 


5 


6.300x10^ 


63 


0.03 


0.01 


0.31 


0.14 


6 


2.160x10^ 


216 


0.03 


0.00 


0.00 


0.36 


7 


7.601x102 


760 


0.04 


0.02 


0.43 


-0.27 


8 


2.724x10^ 


2 725 


0.04 


0.03 


0.60 


0.08 


9 


9.903x10^ 


9 910 


0.05 


0.07 


1.48 


-0.14 


10 


3.644x10^ 


36 446 


0.05 


0.01 


0.21 


0.10 


11 


1.352x105 


135 268 


0.06 


0.04 


0.69 


0.09 


12 


5.056x10^ 


505 861 


0.07 


0.04 


0.66 


-0.04 


13 


1.903x10^ 


1 903 890 


0.08 


0.04 


0.51 


-0.24 


14 


7.205x10^ 


7 204 874 


0.09 


0.01 


0.06 


-0.13 


15 


2.741x10^ 


27 394 666 


0.09 


0.05 


0.49 


-0.33 


16 


1.046x10^ 


104 592 937 


0.09 


0.01 


0.07 


-0.09 


17 


4.009x10^ 


400 795 844 


0.11 


0.03 


0.29 


0.74 


18 


1.543x10^ 


1 540 820 542 


0.12 


0.13 


1.09 


0.44 


19 


5.942x10^ 


5 940 738 676 


0.10 


0.01 


0.15 


0.26 


20 


2.298x10^° 




0.13 






-0.42 


21 


8.895 xlOi° 




0.15 






-0.02 


22 


3.451x10^1 




0.17 






0.62 


23 


1.341x10^2 




0.18 






0.61 


24 


5.228x10^2 




0.20 






1.61 


25 


2.039x10^''^ 




0.19 






-0.04 


26 


7.970x10^3 




0.26 






-0.05 


27 


3.122x10^4 




0.25 






0.00 


28 


1.225x10^5 




0.24 






0.33 


29 


4.831x10^5 




0.28 






0.20 


30 


1.883x10^6 




0.30 






-0.13 



Table 2: Table continued on next page 



17 



N 


Rosenbluth 
estimate 


Exact 
value 


^est 

% error 


True 
%error 


X 




31 


7.426 xlO^^ 




0.33 






0.97 


32 


2.945x10^^ 




0.45 






0.59 


33 


1.160x10^^ 




0.34 






0.19 


34 


4.561x10^^ 




0.47 






0.44 


35 


1.800x10^9 




0.40 






0.23 


36 


7.121x10^9 




0.52 






0.29 


37 


2.823x10^0 




0.57 






0.67 


38 


1.122x10^1 




0.67 






-0.03 


39 


4.417x10^1 




0.65 






0.71 


40 


1.763x10^2 




0.83 






1.30 


41 


6.979x10^2 




0.84 






1.02 


42 


2.738x1023 




0.78 






0.37 


43 


1.088x102^ 




0.82 






-0.16 


44 


4.341x1024 




0.93 






2.12 


45 


1.704x1025 




0.97 






0.52 


46 


6.802x1025 




1.10 






0.73 


47 


2.673x1026 




1.07 






0.41 


48 


1.058x102^ 




1.02 






0.60 


49 


4.209x1027 




1.14 






0.26 


50 


1.664x1028 




1.28 






0.29 



Table 2: Continued:- Degenerate Rosenbluth estimate of the number of lat- 
tice animals of size on a two dimensional square lattice using 2.5 x 10^ 
sample clusters, each grown to = 50 with p = 0.60; exact results from 0; 
calculation of error estimate described in text; 'true' error is fractional dif- 
ference of Rosenbluth estimate and true value; x ^ ^-^e defined in the 
text. 
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